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C^! Abstract 

In this work, the quark hadron phase transition in a chameleon Brans-Dicke model of brane 



world cosmology within an effective model of QCD is investigated. Whereas in CBD model of 



brane world cosmology the Friedmann equation and conservation of density energy are modified, 

X. 

resulting in an increased expansion at early universe. These have important effect on quark hadron 



phase transition. We investigate the evolution of the physical quantities relevant to quantitative 
description of the early times, namely, the energy density, p, temperature, T, and the scale factor, 
o, before, during and after the phase transition. Our analyzes shows that, the quark hadron phase 
transition has took place at approximately nano second after the Big Bang. 
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I. INTRODUCTIONS 



According to the standard model of cosmology, when the early universe expanded and 
cooled, it past through a series of symmetry breaking phase transitions, which predict some 
topological defects. At the early times, above the phase transition, the universe is filled by 
a quark-gluon plasma, (QGP), so study this quark-gluon plasma to hadron phase transition 
will help to provide better understanding of the evolution of the early universe. This early 
universe phase transition has been studied in detailed for over three decades Ufl. Most 
studies have shown that one may obtain first, second, or higher order transition. In addition, 
the possibility of no phase transition was predicted in ?]]. It is well known that the order 
of the phase transition strongly depends on the mass and flavor of the quarks. An early 
study of first order quark-hadron phase transition in the expanding universe is done in [8]. 
As the deconfined quark-gluon plasma cools below the critical temperature Tc =150 MeV, 
it becomes energetically favorable to form color confined hadrons (mainly pions and a tiny 
amount of neutrons and protons, due to the conserved net baryon number). However, this 
new phase does not form immediately. As is characteristic for a first order phase transition, 
some supercooling is needed to overcome the energy expense of forming the surface of the 
bubble and the new hadron phase. When a hadron bubble is nucleated, latent heat is re- 
leased, and a spherical shock wave expands into the surrounding supercooled quark-gluon 
plasma. This reheats the plasma to the critical temperature, preventing further nucleation 
in a region passed through one or more shock fronts. Generally, bubble growth is described 
by deflagrations, with a shock front preceding the actual transition front. The nucleation 
stops when the whole universe has reheated to T c . After that, the hadron bubbles grow 
at the expense of the quark phase and eventually percolate or coalesce. The transition 
ends when all quark-gluon plasma has been converted to hadrons, neglecting possible quark 
nugget production. The physics of the quark-hadron phase transition, as well as the cos- 
mological implications of this process have been extensively discussed in the framework of 
general relativistic cosmology in refs. 



The scalar-tensor theory was conceived originally by Jordan, who embedded a four- 
dimensional curved manifold in five dimensional flat space-time 20j. The scalar-tensor 
models include a scalar field, <ft, whit non-minimal coupling with the geometry in the gravi- 
tational action, as introduced by Brans and Dicke (BD) 21] . Brans-Dicke model have proved 
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to be usefu 
mology 22 



in producing clues to the solutions for some of the outstanding problems in cos- 



231 ] . The mechanism that creates a non- minimal scalar field coupling to the 



geometry can also lead to a coupling between the scalar and matter field. Two such exam- 



24 



251 ] . The scalar field in quintessence 



pies are the quintessence and chameleon scalar field 
mechanism has a very small mass and couples to matter with gravitational strength. The 
authors of ref. 26] introduced the chameleon mechanism in which the scalar field directly 



couples to matter with order unity strength. In this mechanism the mass of scalar field de- 
pends on the local mass density. Indeed, the chameleon proposal provides a way to generat 
an effective mass for a light scalar field via the field self interaction, and the interaction 
between matter and scalar fields. When the chameleon coupling is used in the Brans-Dicke 



model, this called the chameleon-Brans-Dicke model [27J. Recently, the chameleon-Brans- 
Dicke model has been used to describe, holographic dark energy (HDE), agegraphic dark 
energy (ADE), and new agegraphic dark energy (NADE) 281 ] . These investigations have 
concentrate on the large-scale aspects of these models. Moffat et al., [29], have studied 
3 PN-parameter in Jordan-Brans-Dicke cosmology, and the recent work of Perivolaropoulos 
30], and ref.jsjj have considered observational solar system constraints of Brans-Dicke and 
Chameleon-Brans-Dicke models, respectively. 

Over the past decade the possibility that our 4D universe is a brane embedded in a higher 



dimensional space-time has attracted considerable interest [32] . This scenario has been in- 
vestigated for the case in which the bulk is five dimensional and it has been shown that it 
can result in a theory of gravity which mimics purely four-dimensional gravity, both with 
resp ect to the classical gravitational potential and with respect to gravitational radiation 
33j . The existence of branes and the requirement that matter fields should be localized on 
the brane leads to a non-conventional cosmology which has seen a large number of studies. 
Of interest in the present study are the brane-world models in the context of Chameleon 
Brans-Dicke (CBD) gravity. Interestingly, it will be show that in such a scenario and in 
the presence of a CBD field in the bulk, due to nonminimal coupling between scalar field 
and matter field, the conservation equation for the matter field can not be satisfied. It is 
therefore of interest to study the phase transition of quark-gluon plasma to hadrons in the 
context of the CBD brane world theory of gravity. The quark-hadron phase transition in the 
context of conventional brane-world gravity and in Brans-Dicke brane-world gravity have 
been studied in 34 
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II. GENERAL FRAMEWORK AND FIELD EQUATION IN CHAMELEON- 
BRANS-DICKE(CBD) BRANE SCENARIO 



We consider a chameleon-Brans-Dicke action as 

S = -^2 J d 5 xV=g{4>R- < ^d A <f)d A ( j>-V( ( j>)} + J d 5 xf(<P)L m , (1) 

where g is the metric determinant, R is the Ricci scalar constructed from the metric Qab-, <P 
is a CBD scalar field, a; is a dimensionless coupling constant which determines the coupling 
between gravity and the CBD scalar field, L m represents the density Lagrangian for the 
matter fields and V{4>) is the potential. Latin indices denote 5-dimensional components 
(A,B = 1, . . . , 5) and for convenience we choose kL\ = 8nG^ = 1. The last term on the 
right hand said of (1), namely f((f))L m , indicates non- minimal coupling between scalar field 
and matter, where f(<f>) is an arbitrary function of <p. One can obtain the gravitational field 
equation by taking variation of the action (1) with respect to the metric qab and 4> as 

1 „ 1 



5 G AB = Rab - ~9abR = - (rt B + I^Tab) , (2) 



where 



Tab = ^{ Va0V b - ^g AB V 'c<f>V c V>} + {v A V B - SabV^V^} - l -g AB V{(t>), (3) 
and 

W = (3^4) (><*> " 5™)*) + W+T) (I^W " 5™>) ■ « 
Here Tab indicates the energy-momentum tensor of the matter content of the five dimen- 
sional space time and is the trace of it. It is well seen that by f(<fi) = 1 and V(<f>) = 
the above equation reduces to the results of 36] . We consider a 5-dimensional flat metric of 
the following form 

ds 2 = -n 2 (t, y)dt 2 + a 2 (t, y)jijdx' dx j + b 2 (t, y)dy 2 , (5) 

where i, j = 1, 2, 3. We also assume an orbifold symmetry along the fifth direction y = —y. 
We define the energy-momentum tensor as 

T B = T B | fen -\-T B |f,r, (6) 
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where the subscripts br and bu refer to the corresponding energy-momentum tensors in brane 
and bulk respectively. We assume the barne has tension A and filled with ordinary matter 
and the bulk is empty from ordinary matter. So that, the energy momentum tensors are 



rpA 



b \br = -£-diag(-p,p,p,p,0), 
B \ bu = diag(0,0,0,0,0), 



(7) 
(8) 



where it is assumed that the brane is held at y = 0, and 



P = Pm + A, 
p = p m - X, 



(9) 
(10) 



here the subscript m denotes the matter. There are several suggestion for the value of brane 
tension, A. One is that, suggested by the big bang nucleosynthesis, A > 1 MeV 4 [37]. A 



much stronger bounc 
giving A > 10 8 GeV 4 



for A is due to the null results of submillimeter tests of Newton's law, 



38j | . An astrophysical limit on A which is independent of the Newton's 



law and cosmological limits has been studied in 37|, leading to the value A > 5 x 10 8 MeV 4 . 
We assume that the metric ( 1IVDI) . is continuous but their first derivative with respect to y 
is discontinuous, and then their second derivative with respect to y includes the Dirac delta 
function. Making use the metric, one can obtain the non-vanishing component of Einstein 
tensor. The (0,0) components of Einstein equation is 



G, 



00 



a I a + M n 2 fa" ^a' fa' b' 
a \ a b } b 2 V a a V a b 




i (T* + /(0)T OO ) , 



(11) 



where 



T, 



oo 



d b _ u(j) \ fn^ 2 
\ + b 26 + U 



h" + d>' | 3: 



b 



2' 



(12) 



The (i, j) component of Einstein tensor is as 



where 



T; 
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ay 
n) 



ay 
b) 



4> + <j>( 



ou 6 h n 6 N 
— — + 2 h - 

2 a n by 



+ 



n a' n' b' 

+ 2 + - 

a n b 



1 



aV(0) U 



(14) 



The (0,5) and (5,5) components of Einstein tensor are as 



. n' a a' b a' 
G 5 = 3 | + -7 - - 

n a a b a 



l T <t> 

j. 05' 



a' ( a! 



n 



G 55 = 3-{ -- - + - + 



where 



a \ a n 



n 



a I a n 



rr \a \a n 



\ (rl + /(0)t 55 ) , 



(15) 
(16) 



= <t>'~<t> --^ 



n 



T£ 5 = + 3 



d 


h 






a 




n 


20J 





0' 3- + - - 
a n 2 



(17) 
(18) 



And equation of motion of CBD scalar field is 

Ll (n' a' b' 



+ 



d b h 

3- + T-- 

a b n 



n\ 2 
b 



+ 



+ 3-- , 

nab 



+ 



(SH-^ 

(-n 2 ) (|0\/(0) - |\/(0) 
(3a; + 4) 



£19) 



Here, X : = X' := ^ and X := Since, the second derivative of metric consists the 
Dirac delta function, according to [47], one can define 

W" = W" + [W']5(y) (20) 

where W" is the non-distributional part of the double derivative of W(t,y), and \W'\ is the 
jump in the first derivative across y — 0, which defined by 

[W'\ = W'(0 + ) - w\o-). 

The junction functions can be obtained by matching the Dirac delta function in the com- 
ponent of Einstein tensor with the component of brane and scalar field energy-momentum 
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tensor. From the (0,0) and component of field equation we have, respectively 

(p + (u + l)p)/(0 o ) - TP0o/~(0o) ], (21) 

0O/(0O 





-1 


a 6 


(3w + 4)0o 


K] 


1 




(3w + 4)0o 




(3p - p) 


00^0 


(3u> + 4)0o 



(3(w + l)p + (2w + 3)p) /(0o) + (3p - p) | , (22) 

/(0o)-^o/~(0o)|, (23) 



where 



27 = 3w m - I. 

I — I 

For /(0) = 1, these equations reduce to the junction relations in [36|, and the first two 



junction relation are equivalent to [40J|. Using (0, 0) component of the Einstein field equations 
for the brane located at y — and the equations which represent the jump conditions ( l2Ti) . 
(j22p and fl23|) one can find the Friedmann equation as 



/:/ 2 ■ XI// - ^T) = iio _ 1 , { |(3p-p) 2 4 (2 + u.>-)p--u.>pp-2/; 



/ 2 (</>o) 



6 ~ y 4(3w + 4) 2 ^ 2 

+ (3 P - P) 2 (^ - + ^(3p - p) + 2p) (3p - p)0o/(0o)/(0o) 

+ ^(3a; + 4) 2 0on0o)|, (24) 

where if is the Hubbel parameter, T = ~ and the subscript denote the quantity is on 
the brane. From (0,5) component of field equation and using (I24p . (T25]) . fT26|) . we obtain 
continuity equation on the brane as 

p + 3tf (p + p) = -(3p + p)4^</>o (25) 

as we expected the energy, due to the interaction between matter and scalar field, is not 
conserved. Note that in all of relation in this section we select no = 1 and b Q — 1, without 
any less of generality. By substituting the energy density and pressure, the generalized 
continuity equation explicitly confirms the junction condition (|2~Tj) - (I2"3"|) . Using (20) and 
(23) we can obtain the equation of motion for on the brane as 

f Q+ SHT = J^{^(0o)-^f(0o) + (l + f )0o/(0o)/(0o) 

|0O^(0 O ) - |V(0 O ) (2g) 



(3w + 4)0o 



for obtaining (26) we have assumed 0" = 



III. QUARK- HADRON PHASE TRANSITION 

As mentioned in the Introduction, the subject of phase transition from quark-gluon phase 
to hadron phase is a notion fundamental to study of particle physics, particularly in the 
context of lattice gauge theories. Any way, it is an usual part of any study dealing with the 
underlying mechanisms responsible for the evolving universe at its early stages of formation 
in which a soup of quarks and gluons interact and undergo a phase transition to form 
hadrons. Therefore it is important which we have an overview of the basic ideas before 
attempting to use the results obtained from such a phase transition and apply them to the 
study of the evolution of the early universe within the context of the CBD brane-world 
scenario. For further review see (34J and references there in. The equation of state of matter 



in the quark-gluon phase is as [34 ] 



p q = 3a q T 4 + V(T), p q = a q T A -V(T), (27) 



here the subscript q denotes that the matter is quark, a 



g q = 16 + (21/2)N f + 14.25 = 51.25 and V(T) is given as [39] 



tt 2 /90)^, which for N f = 2, 



V(T) = B + lT T 2 - a T T\ (28) 

where B is the bag pressure constant, «t = 77r 2 /20, and jt = f^s/4, which m s is the mass of 
the strange quark in the range m s G (60 — 200) MeV. The form of the potential V corresponds 
to a physical model in which the quark fields are interacting with a chiral field formed with 
the 7r meson field and a scalar field li| . Results obtained in low energy hadron spectroscopy, 
heavy ion collisions and phenomenological fits of light hadron properties give B 1 ^ between 
100 and 200 MeV. In the hadron phase one takes the cosmological fluid to consist of an ideal 
gas of massless pions and of nucleons described by the Maxwell-Boltzmann statistics, with 
energy density p h and pressure ph, respectively. Hence the equation of state is given as 

Ph = T^ph = a n T 4 , (29) 



S 



where a n = (ir 2 /90)gh and gh = 17.25. The critical temperature T c is defined by the condition 
p g (Tc) = ph(Tc) [TJ, and for m s = B l l A = 200 MeV is given by 



I It + \H% + 4:B(a a + «t- a v ) , , 

T c = 4 / rT V 7 V 9 7 — « 125MeV. 30 

y 2(a g + a T - 

Since the phase transition is first order, all the physical quantities, like the energy density, 
pressure and entropy exhibit discontinuities across the critical curve. 



A. Behavior of temperature in CBD brane gravity model during QHPT 

In this section, we will study the quark-hadron phase transition, in chameleon Brans-Dicke 
brane-world scenario. The quantities which we want to trace through the quark-hadron 
phase transition in the chameleon Brans-Dicke brane-world scenario are the temperature T 
and the scale factor a. For obtaining these parameters we use the results which is obtained 
in sec. 2. For simplicity, we take below assumptions 

0o = Na n , (31) 

where N and n are constant. In fact there is no compelling reason for this choice. However, it 



411 ]. It is quite well known now 



has been shown that for small n it leads to consistent results 
that the solar system experiments give the result for the value of u is u > 40000 
Also, this value has obtained in a theoretically study in CBD model of g ravity 



42 



43]. 



m- A 

case of particular interest is that the product \n\u results of order unit y 14 1] , and this show 



431 ]. found that u is 



that the value of n has to be very small~ 10~ 4 . Also the authors of 
smaller than 40000 on the cosmological scales, and they have shown that the evolution of 
the effective gravitational constant is derived by the evolution of the JBD scalar field. They 
have obtained that the constraint of the type u > 500, 1000 would correspond to a variation 
eff ' — — < 1.38% and 0.69% respectively. Also, there is another useful research about this 



Go 



context which is done in 44]]. They obtained the observational constraints on BD model in 
a flat universe with cosmological constant and cold DM using the latest WMAP and SDSS 
data. They found that within 2a range, the value of u satisfies u < —120.0 or u > 97.8. 
They also obtained the constraint on the rate of change of G at present 

G 

- 5.53 x lO^V 1 < — < 3.32 x lO^V 1 (32) 

G 



at 2cr confidence level. So in our case with assumption (1201) we get 



5.53 x 10 



-20„-l 



-20„-l 



< - = ^ = £if < 3.32 x l(r /u s 
G © 



(33) 



This relation can be used to put an upper and lower bound on £. Assuming the present 
value of the Hubble parameter to be Hq — 2.11 x lO -18 ^" 1 , we obtain 

- 0.026 < n < 0.015. (34) 

It is well seen that for upper bound of n the quantity un ~ 200,5, 2.5 for u = 40000, 1000, 500 
respectively. These results have an excellent agreement with 4ll444l]. Therefore, the conser- 
vation equation takes the form below 

p g + 3#(p g +p g + (3p, + p 9 -2A) -/(&,)) =0 (35) 

We assume the CBD brane world is in the quark-gluon phase, so that T > T c . Therefore, 
using ([27D, 02B and (ESD we have 

[il2a q -Aa T )T^ + 2 lT T^f 



H 



1 



(36) 



3 4a g T 4 + ((6a, + 2« T )T 4 - 2 7t T 2 - 2B - 2Aj f /(0 O ) 
and using (27), (28) and (31) the Friedmann equation, (24), takes the following form 

H 2 = —„ 7~7777~~ ~ ~ — 777 X (37) 



4(3w + 4) 2 (l + n- fn 2 )0 2 



-Hf + (2 + 3co + co 2 )(p q + A) - u(p q + X)(p q - A) - 2( Pq - X) 2 / 2 (0 O ) 



where 



+ 



2(p, - A) - l -uH^ ^00/(00)7(00) + ^(3w + 4) 2 0o^(0 o ) 



#1 = (3p,-p,-4A). 

So one can obtain an expression for describing the evolution of temperature of the CBD 
brane universe in the quark phase as 



-3 



T 



4a q T A + (6a, + 2a T )T 4 - 2 7t T 2 - 2B - 2A |/(0 



2(3w + 4)(1 + n - f n 2 )0 2 ((12a, - Aa T )T 3 + 2 lT T 



(38) 





~0J 


•f 






.6 




— UJ 


+ 




~2~ 



Ti + 2T 2 



Ti0o/(0o)/(0o) + 



3a; 
~2~ 

10 



1 T r 



^/ 2 (0o) 



+ -(3 W + 4) 2 o ^(0o) 



where 



Ti = (4ct T T 4 - 4 7T T 2 - AB - AX), 
T 2 = ({a q + a T )T A - lT T 2 -B-X 
T 3 = ((a, - a T )T A + lT T 2 + B + X s j. 



B. Behavior of temperature in CBD brane gravity during QHPT with V(T) = B 



In this case the equation of state of the quark matter is given by the bag model, namely 
Pq — (Pg ~ 4B)/3. Actually, the formula of this subsection are as same as previous section 
with 7t = «y = 0. So the relation (33) becomes 



H 



12a q T 3 T 



12a q T A + 6 ( 3a q T A - B - X ) nf((f) ) 



(39) 



From this relation we should obtain scalar factor with respect to T. Also f l3Tj) takes the 
following form 



T 



Aa q T 4 + 2 (3<T 4 - B - a) |/(0 o ) 



la q {3u + 4)(1 + n - |^ 2 )^T 3 



(40) 



x 



+ 



^T 2 + (2 + 3w + w 2 )T 3 2 - c^3<T 4 + 5 + AjT 2 - 2T 2 2 J / 2 (0 O ) 



— a; 

Ti + 2T 2 

2 



where 



T 1 = -4(B + A), 



T 2 
T 3 



a q T 4 -B-X 



a q T 4 + B + X 



In the above relation scalar field is a function of scalar factor and scalar factor is a function 
of temperature T which is obtained from relation ( 1391) . 
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C. Growing up the volume fraction of hadron in CBD brane gravity during QHPT 



During the quark-hadron phase transition, p q (t) decreases from p q (T c ) = pq to Ph(T c ) = 
Ph but, the temperature and pressure are constant. Also the quantities like the enthalpy 
W — (p q + p q )a 3 and entropy S = sa 3 are conserved. In phase transition temperature, 



T c = 125 MeV we have p q « 5 x 10 9 MeV 4 , p h « 1.38 x 10 9 MeV 4 and 



in which /v is constant during the phase transition. Following {?], |31. 



Pc 



« 4.6 x 10 8 MeV 4 
, one can replace 



p(t) by h(t), the volume fraction of matter in the hadron phase, by defining 

p q {t) = (l + mh{T))p Q , 



(41) 



here m = puj Pq — 1 = constant. At the beginning of the phase transition h(r c ) = 0, where 
r c is the time corresponding to the beginning of the phase transition, and p(r c ) = Pq, while 
at the end of the transition /i(r^) = 1, where is the time at which the phase transition 
ends, corresponding to p{rh) = pu- For r > the universe enters in the hadronic phase. 
Then with the help of conservation relation we arrive at 



H 



1 



rh 



(42) 



where r is defined as 



Ph 



— = constant. 
Pc + Pq 

The evaluation of hadron fraction during the phase transition is given by equation 



(43) 



l + rh + 





3 




r 


'us 


X < 










.6 






-UJ 


+ 








Y 



3p c +PQ-2\ 
Pc+PQ 



+ rMfiv>o) 



where 



2r(3w + 4)(l + n- %n 2 )(pl 



A x + (2 + 3lu + lu 2 )A 2 2 - ujA 2 (p c - A) - 2(p c - X)' 



/ 2 (0o) 



(44) 



-fA, + 2(p c - A)] AM fa) f fa) + (y)^^ + jj(3c + AfaV{<j>) 



A x = (3p c - pq - (p H - p Q )h - 4A), 
A 2 = (pq + (ph - Po)h + A). 
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D. Behavior of temperature in CBD brane gravity in pure hadronic era 



In this stage of universe, we encounter with a pure hadronic era with following equation 
of state 

Ph = 3ph = 3a n T 4 , (45) 
then with the help of conservation relation one can obtain 

1 12a 7r T 3 T 



H 3 4a.T4+(2p£-2A)f/(0) : 
and from the Friedmann equation one arrive at 



(46) 



T 



-3 ^4a,T 4 +(2p b h -2A)f f(cf>) 
18a^T 3 (3w+4)(l+n-|n 2 )0g 



+ (2 + 3co + co 2 )(3a n T A + A) 2 - u(3a n T 4 + X)(a w T A - A) 



-2(a w T 4 - A) 2 / 2 (0 O ) - 4A - 8a;A 2 + 2(a^T 4 - A) M (<Po) f (fo) 



+16A 2 (|^ - l) + 2 {3u + 4 ) VoV(0 o ) }. 



(47) 



IV. TYPICAL EXAMPLE 

In this section we study our model for some example of /(</>) and V(<p). Since we study the 
chameleon kind of Brans-Dicke theory in brane world scenario, then we need a monotically 



decreasing 



unction for potential where exhibit self interaction and satisfied the chameleon 



costraint 26]. Therefore we consider power law potential as 



n$) = M 5 f^y 



(48) 



where £ > and M is a constant with mass dimension and Ai is dimensionless. In this 
case we assume that the scalar field is dimensionless. The authors of 26j have found that for 
small value of £ G (0, 2) the magnitude of M and Ai are ~ 10 -3 eV. Therefore, with point 



of view the above mention and use of (31), the potential takes the following form 

10 -(3£+12) 10 -(3f+12) 



(49) 



Here a is the scale factor of the universe. To obtain an explicit view of the universe we have 
to find the scale factor. Since other quantities in this theory are versus MeV, then the scalar 
potential part is very small with respect to other terms and we can withdraw it. Also solving 
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the basic equation of our model for an arbitrary function of f(<f>) is too complicated, so that 
we complete our investigation with f(<f>) = <fi and we will obtain a and T for all subsection 
in the pervious section. 



A. Bhavior temperature in CBD brane world during QHPT 

For obtaining the relevant quantity we use the results which is obtained in sec. 2. So that, 
setting f((f>) = (f> and /((/>) = 1 in equations (32), we find the conservation equation as 

p + 3H^p q + Pq + (3 Pq + p q -2X) 7 ^j = 0, (50) 

and the Hubble parameter takes the form 

2 ((6a, - 2a T )T 3 + 1t t)t 



H= a - 

a 



{a q (3n + 12) + na T )T 4 - n-f T T 2 - n(B + A) ' 
2(aT 3 + 7T T)T 



&T 4 - n lT T 2 - B 1 



(51) 



where 



a = Qa q — 2dT, 
b = (12 + 3n)a q + nar, 
B 1 = n(B + X), 



integrating the above equation gives 



a(T) = C^bT 4 - n lT T 2 - B - X 



(6+3n)ag + ri 



exp 



7 



2^/bB~ Q 



arctan 



-DO 



(52) 
(53) 
(54) 



»7r 
2b 



where C is a constant of integration and the other constants are as follows 



Bo = B 1 + 



n 2 7 2 
~46~' 



na 

7 = ( ^ + 1)7t- 
Also we obtain the Friedmann equation as follows 
1 f/w-2 



4(3w + 4) 2 # 



-fT 1 + 2(p,-A) 



(55) 

(56) 
(57) 

(58) 



+ 



u> 



Hi + (2 + 3cu + w 2 )(p 9 + A) 2 - w(p 9 + A)(p, - A) - 2{ Pq - Xf 
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where 



#i = (3p 9 -p 9 -4A), 



9 



l + n 



So one can obtain an expression for describing the evolution of temperature of the CBD 
brane universe in the quark phase by combining the above equations as 

(6T 4 -n 7T T 2 - B x ) ffuJ-2 



2(3w + 4)(aT 3 + 7 r T)v^\^ 4 



-H 1 + 2(p q -X))H 1 



+ 



to 



-Hi + (2 + 3oo + co 2 )(p q + A) - co(p q + X)(p q - A) - 2(p q - \)< 



(59) 



Substituting (l2Tj) and (|28|) into (1591) and using (}3~Tj) and (|55|) . we obtain T as a function of 
r. We solved it numerically and the result is plotted in FIG.l. Fig. la shows the behavior of 
temperature as a function of the cosmic time, r, in a CBD brane-world gravity, during the 
quark hadron phase transition for different values of u> with n = 0.015 , Nq = 2 x 10 5 and 
A = 10 x 10 8 MeV 4 . This plot shows that the phase transition accrued at approximately 
nano second after Bing Bang. FIG. lb shows the behavior of scalar factor as a function of 
the temperature, T, in a CBD brane-world model of gravity, during the quark hadron phase 
transition. This figure shows, when phase transition took place, the effective temperature 
of the universe is decreased and scale factor is increased. This means that our model satisfy 
the expansion of universe during the quark hadron phase transition at approximately 10 -9 
s after the Bing Bang. 

B. Behavior of temperature in CBD brane gravity during QHPH with V(T) = B 



The formula of this subsection are as same as previous section with 
setting / = in (36) we have 



H 



12a q T 3 T 



3(4 + n)a q T 4 - (B + X)n' 



o>t = 0. So 



(60) 



Integrating (57) gives 



a(T) = c ( 3(n + 4)T 4 - n{B + A) 



4 

"3(n+4) 



(61) 
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(a) (b) 

FIG. 1: (a): shows the temperature of the universe during the quark hadron phase transition for 
different values of u = 10 4 (solid), 1.5 x 10 4 (dashed), 2 x 10 4 (dotted), 2.5 x 10 4 (dotted - dashed) 
versus cosmic time(b): shows the scale factor of the universe versus temperature of the universe in 
quark gluon plasma phase with V(T) ^ constant versus temperature. We have taken N = 2x 10 5 , 
A = 10 9 , n = 0.015 and B 1 / 4 = 200 Mev. 



and 



T 



3(4 + n)a q T 4 - (B + X)n 



+ 



2Aa q T 3 (3uj + 4)V26 
8u(B + A) 2 



x 



4 ( (2 - u){B + A) + 2{p q - A) ) {B + A) 



(2 + 3uj + io 2 )(p q + A) - u(p q + A)(p, - A) - 2( Pq - \)' 



Using (!3T|) and (IUT|) and substitute (1271) . (!28|) into f|62|) . one can obtain the quantity T as 
a function of T and r. We solved (59) numerically and the final result is plotted in FIG. 2. 
which shows the behavior of temperature as a function of the cosmic time r in a CBD 
brane-world filled with quark matter for different values of u, for V(T) = B. This figure 
confirm the result of FIG.l 



C. Growing up the volume fraction of hadron in CBD brane gravity during the 
QHPT 

In this subsection we want to consider the obtained results in C subsection in the pervious 
section for /(0) = <p an d £ = 1- As was mentioned ago the pressure during the phase 
transition is constant, p c ~ 4.6 x 10 8 MeV 4 and also density of quark matter and hadron 
matter at the transition time are p q ~ 5 x 10 9 MeV 4 , ph ~ 1.38 x 10 9 MeV 4 respectively. 
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50 
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100 120 140 160 ISO 200 

T T 

(a) (b) 

FIG. 2: (a): shows the temperature of the universe during the quark hadron phase transition for 
different values of oj = 8 x 10 3 (solid), 1.4 x 10 4 (dashed), 2 x 10 4 (dotted), 2.6 x 10 4 (dotted - dashed) 
versus cosmic time for V(T) = B. (b): shows the scale factor of the universe versus temperature 
of the universe in quark gluon plasma phase versus temperature. We have taken N = 2 x 10 5 , 
A = 10 9 , n = 0.015 and V(T) = £ 1/4 = 200 Mev 



Therefore, in this case (39) reduces to 

1 



H 



rh 



where r is denned as 



3 1 + r/l+ ( 3pc+PQ-2X j\n 
V Pc+PQ J 2 



r = P -^Q = -5.6 x ICT 6 . 



(62) 



(63) 



Pc + PQ 

Integrating of (60) one can obtain relation between the scale factor on the brane and the 
hadronic fraction h(r) as 



a(^a(r c )(l^(l + 2 ^ + (l + ^)) 



, (64) 

pc + p q ; ■ x r ! 

here we assume h(r c ) = 0. So, using equations fl58|) we obtain the time evolution of the 
matter fraction in the hadronic phase as 



h = -3 x 



(^h 1 +2(j, c -\fjh 1 



2r(3cu + 4)V^ 



+ 



L6 



hi + (2 + 3a; + u 2 ){ Pq + A) 2 - cu(p q + A)(p c - A) - 2(p c - Xf 



where 



h = (3p c - p q - 4A), 
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Pq = PQ(l + mh(T)). 

After substituting f[3"Tj) and (JMI) into (fB"5]) . the hadron fraction, /i(t), as a function of r is 
obtained. The numerical results of h(r) is presented in FIG. 3 for different value of u. FIG. 3a 
shows the volume fraction of hadron during the quark hadron phase transition for different 
values of u = 10 4 (solid), 1.5 x 10 4 (dashed), 2 x 10 4 (dotted), 2.5 x 10 4 (dotted - dashed) versus 
cosmic time. FIG3.b shows the scale factor of the universe during the QHPT versus volume 
hadron fraction. It is well known that when QHPT has took place the density of quark 
gluon plasma is decreased and the volume fraction of hadron and the scale factor of universe 
are increased. The FIG. 3 shows that this fact has took place at 10 -9 and confirm the results 
of FIG.l too. 




(a) (b) 

FIG. 3: (a): shows the volume fraction of hadron during the quark hadron phase transition for 
different values of u = 10 4 (solid), 1.5 x 10 4 (dashed), 2 x 10 4 (dotted), 2.5 x 10 4 (dotted - dashed) 
versus cosmic time, (b): shows the scale factor of the universe during quark hadron phase transition 
versus volume fraction. We have taken N = 20000, A = 10 9 , n = 0.015 and B 1 / 4 = 200 Mev 



D. Behavior of temperature in CBD brane gravity after QHPT 

In this stage of universe we encounter with a pure hadronic era with following equation 
of state 

Ph = 3ph = 3a w T 4 , (65) 
then with the help of conservation relation one can obtain 

1 Yla^f 
~3 4^T 4 + (2p£-2A)f (66) 
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integrating of (|66|) leads to 

a(T) = a (T c ) (4 + ^r^a^T 4 — \n 
and from the Friedmann equation one arrive at 



i 

" 4+3n 



(67) 



T 



({A + 3n)a n T 4 - n\ 



+ 



4a^T 3 (3w + 4) v / 20 
IX 2 



4A(wA - 2a w T 4 ) 



(68) 



+ (2 + 3oo + w 2 )(3a w T 4 + A) - w(3a w T 4 + A)(a w T 4 - A) - 2(a^T 4 - A)' 



FIG. 4 shows variation of the temperature of hadronic fluid filled CBD brane- world as a 
function of the cosmic time r for different values of u, with n = 0.015 , Nq = 2 x 10 5 and 
A = 10 x 10 8 MeV 4 . 
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FIG. 4: (a): shows the temperature of the universe after the quark hadron phase transition for 
different values of u = 10 4 (solid), 1.55 x 10 4 (dashed), 2.1 x 10 4 (dotted), 2.55 x 10 3 (dotted - dashed) 
versus cosmic time, (b): shows the scale factor of the universe versus temperature of the universe in 
hadron phase versus temperature. We have taken N = 2 x 10 5 , A = 10 9 , n = 0.03 and B 1 ^ = 200 
Mev 

In FIG. 5 we have plotted the effective temperature and volume fraction of hadron as a 
function of cosmic time. FIG. 5a shows T vs r for different pairs of (u,n) in which um = 1. 
It is seen that all curves are similar to each others and overlap completely. FIG. 5b shows 
h(r) vs r for different pairs of (oJ,n) in which um = 1 and has the same behavior. 

The effective temperature has been plotted in GIG. 6 for different kind of (ui, n) pairs that 
um = c. FIG6.a shows T vs r for (u, ri) = (10 5 , 10~ 3 ), (3 x 10 4 , 10~73) in which un = 100. 
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FIG. 5: (a): shows that the dependence of temperature of the universe during the quark 
hadron phase transition for different values of (cj,n) = (10 4 , 0.0001), (3 x 10 4 , 0.0001/3), (5 x 
10 4 , 0.0001/5), (7 x 10 4 , 0.0001/7) in which oj x n = 1(b): shows the volume fraction of hadron vs r 
for different pairs of(w, n) = (10 4 , 0.0001), (3 x 10 4 , 0.0001/3), (5 x 10 4 , 0.0001/5), (7 x 10 4 , 0.0001/7) 
in which u X n = 1. We have taken N = 2 x 10 5 , A = 10 9 and 5 1 / 4 = 200 Mev 




2.xl0-4 xl0"d°xl0^ xl0-l°xl0- 9 2. x 10 _ 4°x 10 _ d°x 10"^°x 10 _ 1° x 10 -9 



FIG. 6: (a): shows that the dependence of temperature of the universe during the quark hadron 
phase transition for (u),n) = (10 5 , 10 -3 ), (3 * 10 4 , 10" 2 / 3 ) in whicn 

uj x n = 1(b): shows the 

temperature of the universe during the quark hadron phase transition for different values of (oj, n) = 
(10 6 , 10~ 3 )(solid), (3 x 10 5 , 0.01/3) (dashed) in which u x n = 1000. We have taken N = 2 x 10 5 , 
A = 10 9 and B 1 / 4 = 200 Mev 

This figure indicates these curves have behavior like FIG. 5a. But FIG6.b shows T vs r for 
(u),n) = (10 6 , 10~ 3 ), (3 x 10 5 , 10 _2 /3) in which un = 1000. It is seen that these curves have 
different behavior. We analyzed our model for another values of c. Our considerations show 
that for n < 0.015, c ^ 120 and 100 ^ u ^ 10 6 , all curves have the same behavior versus 



T 



T 



(a) 
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cosmic time and for c > 120 the temperatures curve have not similar behavior. These results 
are in a good agreement with 41]. 



V. CONCLUSION 

In this paper, we have studied the quark-hadron phase transition in a chameleon Brans- 
Dicke brane cosmological model. We investigated the evolution of the physical quantities 
relevant to the physical description of the early times such as; the energy density, temper- 
ature and scale factor, before, during, and after the phase transition. We found that for 
different values of n < 0.015 and 100 ^ u ^ 10 6 phase transition occur and with increas- 
ing time the effective temperature of the quark-gluon plasma and the hadronic fluid will 
be decrease. We plot the effective temperature and scale factor of universe on the differ- 
ent stage of phase transition for various values of u and n. All of our plots, show that 
the effective temperature and the FLRW scale factor decrease and increase respectively by 



possing time. We then compared our results with the results presented in [34M36| . In [34 1 
the authors studied quark hadron phase transition in a RandallSundrum brane model and 
showed that for different values of the barne tension cu and A, phase transition occurs about 



at 10 6 s after the Big Bang. Also in 35J, |36), the authors investigated the quark-hadron 



phase transition in a brane-world scenario where the localization of matter on the brane is 
achieved through the action of a confining potential and showed that for different values of 
the parameters in their model, phase transition takes place. They found that for various 
value of to the phase transition is took place aboutl0 _6 s, but our investigation shows that 
the quark- hadron phase transition has took place at approximately 10 _9 s after Big Bang. 
This is a difference between the results of our study and the studies of other researchers. 
This means that due to the interaction between scalar field and matter, modified conserva- 
tion equation of energy-momentum, the relation between temperature and cosmic time is 
changed, then the rate of expansion of the universe is increased at the early times. We plot 
the effective temperature versus cosmic time for cases which nu = c. These figures show 
that for n ^ 0.015 and u G (10 2 , 10 6 ) in which un ^ 120 all curves completely have similar 
behavior versus cosmic time, and by increasing the c, decreasing the effective temperature 
is faster and the curves with different pairs of (u, n) have not similar behavior versus cosmic 
time. Indeed, this means that we have to be urn ph 1 in the universe. 
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